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Abstract. We present the 0.1-200 keV BeppoSAX spectrum of MR 2251-178 observed at two epochs in 1998 
separated by 5 months. Both epochs show identical spectral shape and X-ray flux. Analysis of the combined 
spectra allow us to confirm the presence of the ionized Fe Kq line detected by ASCA and to test the presence 
of reflection from ionized material. The best self-consistent spectral fit is obtained when including a contribution 
from a mildly ionized reflector (^o.oi-ioo ~ 1625 erg cm s~^, i.e. Uo.i-io ~ 0.14) with a reflection fraction 
Rion ~ 0.11. An exponential cut-off to the direct power-law continuum is then required at E ~ 100 keV. The low 
energy spectrum is absorbed by ionized matter with a column density Nw ~ (S.l'Ji^'}^) x 10^^ cm^^ and an X-ray 
ionization parameter Uo.i-io ~ 0.06 ±0.01. The warm absorber column is slightly lower than that measured with 
ASCA. This change could be caused by bulk motion. The BeppoSAX absorber ionization parameter globally agrees 
with the "U versus Flux" correlation found for MR 2251-178 with EXOSAT, Gmga and ASCA. This suggests that 
the absorber is in a state of instantaneous ionization equilibrium. 
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1. Introduction 

The large variety of observed spectral and temporal vari- 
ability properties of AGN is at the origin of the com- 
plex AGN classification scheme. Attempts to find a sim- 
ple explanation for the different classes have resulted in 
the "unified model" of AGN (Antonucci [1993D , in which 
the viewing angle of the observer is the key parameter. 
AGN luminosity is an additional parameter which can di- 
rectly influence the physical state of matter located in 
the innermost regions of the nucleus. However the role of 
AGN luminosity is at present still not well understood. 
For instance, recent X-ray studies are still inconclusive 
about the existence of any fundamental differences be- 
tween Radio Quiet (RQ) QSOs and RQ Seyfert 1 galax- 
ies (Sy 1). These two classes have often been regarded as 
high- and low-luminosity manifestations of a same phe- 
nomenon (e.g. Staubert & Maisack 1996; Lawson & Turner 



properties appear to be comparable in RQ QSOs and Sy 
Is (George et al. ibid.). However, the strong Compton re- 
flection "hump" which is common in Sy Is has not been 
clearly detected in luminous RQ QSOs (Williams et al. 
1992| , Laws on & Turner |1997|, Ree ves et al. |1997| , George 
et al. |2000| , Reeves & Turner pOOOJ ). An analy sis of ASCA 



1997). For example, 2-10 keV power-law spectral slopes of 



\ow-z RQ QSOs have been found to be similar to th ose of 
low luminosity Sy Is (F ^ 1.9 — 2, e.g. George et al. 200C , 
Reeves & Turner [2000 ). Variability and warm absorber 
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data for high z RQ QSOs (Vignali et al. 1999) shows that 
the 2-30 keV (source frame) spectra are well described by 
a single power-law. Nandra et al. ( 1997b| ) showed that Fe 
Ka line emission in RQ QSOs is weaker than in Sy Is and 
that the line profile and centroid energy both seem to de- 
pend on the AGN luminosity. This trend has recently been 
confirmed by George et al. ( 200C| ). Several mechanisms 
have been proposed to explain the differences between 
RQ QSOs and Sy Is. For instance relativistic beaming 
could enhance the underlying X-ray continuum emission 
but not the iron line, which is thought to arise by fluores- 
cence from cold or hot reflecting material (Williams et al. 
1992 ). Another interpretation is that the X-ray luminos- 
ity of the central source regulates the ionization state of 
the Compton scattering medium which is probably the ac- 
cretion disk itself. As the luminosity increases the atoms 
in the disk become increasingly ionized. At high ioniza- 
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tion states the iron line is emitted a t highe r energies and 
its flux may d ecreas e (Nandra et al. 1997b and r eference s 
therein; Matt |l998t Ross et al. 1999, hereafter |RFY99|) . 
Ultimately, the optical depth for photo-electric absorp- 
tion decreases to such an extent that the continuum flux 
<30 keV is mostly reflected from the ioniz ed disk and no 
longer absorbed (Basko et al. 1974 RFY99). Narrow band 
X-ray observatories have failed to detect any "hard tails" 
in RQ QSOs and the reflection component in these sources 
may so far have been confused with the underlying contin- 
uum emission. However, at high X-ray energies the cut-off 
due to Compton recoil and Klein-Nishina effects may be 
detectable. 

In this paper we present broad band BeppoSAX spec- 
troscopy of a RQ QSO. We show that this source has an 
ionized iron line and exhibits a high energy cut-off, despite 
lacking a Compton reflection hump. We compare our data 
with the predictions of reflection models in order to test 
for the presence of an ionized accretion disk. 

MR 2251-178 is a nearby (z = 0.064) X-ray bright RQ 
QSO (Lx ~ 10^^ ergs~^). Early X-ray observations of this 
source with the Einstein Observatory revealed both a soft 
excess and variable intrinsic absorption. This was t he firs t 
detection of a warm absorber in an AGN (Halpern 1984 ) . 
The source was later monitored by EXOSAT between 1983 
and 1984 with a sequence of 15 observations (Pan et al. 



199C ) . During this period the 2-10 keV flux varied between 
1.8 and 3.3 xlO^^^ ergs""'^ cm~^. The changes occurred on 
time-scales as short as ~10 d. In the same period the col- 
umn density of the warm absorber also varied. Strong cor- 
relations between the soft excess flux, the warm absorber 
column density variability and the source luminosity were 
found. A subsequent reanalysis of the EXOSAT data to- 
gether with Ginga observations (Mineo & Stewart 1993) 
showed that the spectral behavior can be well described 
using a self-consistent warm absorber model, and that the 
ionization parameter of the gas is clearly correlated with 
the continuum flux. ASCA observations in 1993 and 1996 



(Otani et al. 1997; Reeves et al. 1997) reveal a weak emis- 



sion line at 6.57_Qgy keV due to partially ionized iron. In 
this case too, the ionization state of the warm absorber 
correlates well with the measured X-ray flux. High ion- 
ization absorpt ion li nes have been observed in the UV by 
Monier et al. ( 200l| ), these appear to be related to the 
X-ray warm absorber. 

The BeppoSAX observations of MR 2251-178 extend 
the spectral range of previous X-ray spectroscopy of this 
object. At the same time they also expand the luminosity 
range of "warm-absorbed" AGN observed by BeppoSAX. 
A preliminary analysis of the warm absorber data is given 
by Orr et al. ( |2000| ). 



2. Observations 

We present results obtained with the Low-Energy 
Concentrator Spectrometer (LEGS; 0.1-4 keV; Parmar et 
al. 1997| ), the Medium-Energy Concentrator Spectrometer 
(MEGS; 1.8-10 keV; Boella et al. |1997|), the High Pressure 



Gas ScintiUation Proportional Counter (HPGSPC; 4- 
120 keV; Manzo et al. |l997| ) and the Phoswich Detection 
System (PDS; 15-300 keV; Frontera et al. |l997| ) on-board 
BeppoSAX. 

All these instruments are co-aligned and collectively 
referred to as the Narrow Field Instruments, or NFL The 
MEGS consists of two grazing incidence telescopes with 
imaging gas scintillation proportional counters in their fo- 
cal planes. The LEGS uses an identical concentrator sys- 
tem as the MEGS, but utilizes an ultra-thin entrance win- 
dow and a driftless configuration to extend the low-energy 
response to 0.1 keV. The non-imaging HPGSPC consists 
of a single unit with a collimator that remained on-source 
during the entire observation. The non-imaging PDS con- 
sists of four independent units arranged in pairs each hav- 
ing a separate collimator. Each collimator was alterna- 
tively rocked on- and off-source during the observation. 

The region of sky containing MR 2251-178 was ob- 
served by BeppoSAX at two epochs: from 1998 June 14 
17:58 UT to June 18 05:33 UT and from 1998 November 
12 21:23 UT to November 16 05:32. Good data were se- 
lected from intervals when the elevation angle above the 
Earth's limb was >4° and when the instrument configu- 
rations were nominal, using the SAXDAS 2.0.0 data anal- 
ysis package. The standard PDS collimator dwell time of 
96 s for each on- and off-source position was used together 
with a rocking angle of 210'. LEGS and MEGS data were 
extracted centered on the position of MR 2251-178 using 
radii of 8' and 4', respectively. The June exposures in 
the LEGS, MEGS, HPGSPC, and PDS instruments are 
60.3 ks, 82.7 ks, 53.3 ks, and 41.1 ks, respectively. The 
November exposures are in the same order 30.5 ks, 61.2 ks, 
13.9 ks, and 35.2 ks. 

Background subtraction for the LEGS and MEGS were 
performed using the standard background files. For the 
HPGSPC it was carried out using data obtained when the 
instrument was observing the dark Earth. Finally the PDS 
background was estimated from the offset field according 
to the standard procedure. 

3. Results 

3.1. Light curves and variability 

We have examined the LEGS 0.1-2 keV and the MEGS 
2-10 keV light curves for the 1998 June and November 
observations as well as the corresponding hardness ratios. 
The data were binned into intervals of 10 ks. The count 
rates and the hardness ratios are both relatively similar 
at the two epochs. The 1998 June light curves are well fit 
with a constant count rate. This gives, respectively, for the 
soft and hard X-ray bands 0.173 ±0.003 and 0.417 ±0.004 
counts s~^ with corresponding fit statistics xt = 1-1 and 
0.9 (degrees of freedom, dof=29). Fits to the November 
data give 0.183 ± 0.004 and 0.433 ± 0.005 counts s-\ re- 
spectively, with xl = 1-5 (dof=24) and xl = 3.4 (dof=21). 
The poor statistics of the constant fits to the November 
data are an indication of variability. In fact, the November 
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Fig. 1. 0.1-2 keV and 2-10 keV light curves and hardness ratios of the June and November 1998 BeppoSAX observa- 
tions of MR 2251-178 



light curves are better fitted with slowly decreasing linear 
functions (minus ^^15% in 70 hours). 

Finally, an interval of 30 ks beginning 10 ks after the 
observation start was ignored in the MECS November 
light curve. This time interval covers a drop in MECS 
count rate due to a change in pointing which could not 
be corrected by using the standard housekeeping data. 
The cleaned and linearized MECS November image shows 
what appears to be an extended, oblong source at the 
position of MR 2251-178. A comparision of light curves 
for different event extraction radii shows that this is an 
artefact caused by the displacement of the source in the 
MECS field of view. A similar displacement occurred in 
the LECS, however the available housekeeping data allow 
this effect to be corrected. 

The 1998 June 2-10 keV flux is (4.03l|^;i^) x 10"" 
erg s"^ cm"2 (MECS) and the November flux is (4.19 ± 
0.42) X 10~" erg s'^ cm'^. 

3.2. Spectral fits and interpretation 

The June and November spectra of MR 2251-178 were 
investigated both individually and as a single, summed 
spectrum by simultaneously fitting data from all the 
BeppoSAX NFI with the help of the XSPEC 11.01 spec- 
tral fitting package. The LECS and MECS spectra were 
rebinned to oversample the full width half maximum 
(FWHM) of the energy resolution by a factor 3 and to 
have, additionally, a minimum of 20 counts per bin to al- 
low use of the x^ statistic. The HPGSPC spectrum was 
rebinned to have a minimum of 40 counts per bin and 
the PDS data were rebinned using the standard logarith- 
mic binning recommended for this instrument. Data were 
selected in the energy ranges 0.1-4.0 keV (LECS), 1.8- 
10 keV (MECS), 8.0-12 keV (HPGSPC), and 15-200 keV 
(PDS) where the instrument responses are well deter- 
mined and sufficient counts obtained. The photo-electric 
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Fig. 2. The summed June and November 1998 BeppoSAX 
spectra of MR 2251-178. The MECS data are represented 
by small filled squares, the HPGSPC data by open circles. 
The fit model is RFY99 with a warm absorber 



absorption cross sections of Morrisson & McCammon 
(1983) and the solar abundances of Anders & Grevesse 



(1989) are used throughout. 

Factors were included in the spectral fitting to allow for 
fiux normalization uncertainties between the instruments. 
For all instruments factors were constrained to be within 



their usual ranges during the fitting (see Fiore et al. 1999). 
In this paper all uncertainties are 90% confidence intervals 
for one parameter of interest, unless stated otherwise. 

Because of the lower exposure times during the 
November 1998 observation the uncertainties on the fitted 
parameters are slightly larger than for those of the June 
spectrum. However, our fit results indicate that absolutely 
no significant spectral change occurred between the two 
epochs which are separated by ~5 months. All the fit pa- 
rameters for the individual spectra are compatible within 
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Table 1. Fits to the summed June and November 1998 spectrum. The fitted energy ranges for the different instruments 
are 0.1-4 keV (LEGS), 1.8-10 keV (MEGS), 8-12 keV (HPGSPG) and 15-200 (PDS). AU models include galactic 
neutral absorption Nh = 2.78 x 10^*^ cm~^ and, apart from the simple power-law (PL), also 3 absorption edges. The 
edge parameters are only listed for the cut-off power-law (COPL) . All energies are given in the source rest-frame. The 
pexriv model is composed of both the reflected and the direct component. In this model the temperature of the 
ionized reflector has been flxed at 10^ K. The parameters without uncertainty intervals were frozen during the flt. 
The ionization parameter is defined as in RFY99, integrating over the energy range 0.01-100 keV. The units of the 
ionization parameter ^ are erg s^^ cm 
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Fig. 3. Gontour plot of the parameters Ecut-off and F in 
the ionized reflection model pexriv. The contours indicate 
68, 90 and 99% confidence levels 
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Fig. 4. Data-to-model ratio for a simple power-law with 
galactic neutral absorption. The MEGS data are shown 
with solid squares and the HPGSPG with open circles. The 
warm absorption and the high energy cut-off are evident 



the statistical uncertainties for any of the models tested. 
For this reason and because the observed fluxes in June 
and November are nearly identical with no large structure 
change in the light curves we have summed the spectra for 
the two epochs in order to obtain a higher signal to noise 
ratio. 

We tried fitting the summed continuum with several 
models (see Table H), starting with a power-law which 
gives a totally unacceptable fit {xi — 4.2, dof=162). 
Strong deviations appear in the ranges E~0.4-3 keV and 
at E^ 80 keV as illustrated in Fig. |[ The spectral com- 
plexity in the soft X-rays is due to the warm absorber and 



is discussed further on. At high energies the data points 
clearly fall below a simple power-law. Highly significant 
improvement is brought with an exponential high energy 
cut-off to the power-law at Ecut-off > 70 keV, (the con- 
fidence level of the improvement is >99.9%; F-statistic 
F>40 for Adof=l, with respect to a simple power-law, 
with or without an Fe Ka line or a warm absorber; the 
value of the cut-off energy is discussed below). However, 
as described in the following sections, the best fit to the 
summed 0.1-200 keV BeppoSAX spectrum of MR 2251- 
178 is obtained with an ionized reflection model. Such a 
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model predicts a drop in high X-ray flux but, unlike the 
case of reflection onto cold matter, a weaker or absent 
"Compton hump" . 

3.2.1. Significance of the high-energy cut-off 

The knowledge of the spectral shape above ^10 keV is 
an important key to the interpretation of the present 
BeppoSAX data. Therefore we have taken special care in 
verifying the quality of the PDS data ( for d etails about 
the PDS instrument see Frontera et al. 



1997) 



It has been shown (by Guainazzi and Matteuzzi, 1997) 
that the total systematics in the PDS 20-200 keV range 
due to background subtraction are lower than 0.026 s~^ 
(corresonding to approximately 5.8 % of the PDS 20-200 
keV source count rate for MR 2251-178, which is 0.45 cts 
s~^). The observed flux of a one mCrab source with pho- 
ton index 1.5-2, is well above the systematics (^ la) up 
to 60 keV and remains above the systematics (> la) up 
to 100 keV). (MR 2251-178 has a flux of ~3 mCrab in the 
20-200 keV range). 

We have also checked the consistency of the spectra 
measured by the two PDS half-arrays (A & B) which see 
the sky (source and background field) through two dif- 
ferent collimators. This was done by fitting either one of 
the PDS half-arrays (A & B) spectra together with the 
LEGS, MEGS and HPGSPG data. Because the data were 
obtained during two epochs (June and November 1998), 
we first tested the half array data for June only and then 
for the summed June and November. This is justified be- 
cause the PDS rocking collimator positions happen to co- 
incide within 0.6 degrees. Using the "cut-off power-law + 
3 edges -|- Fe Ka line" model in Table |l] one finds that the 
cut-off energies are all compatible with those found using 
the combined PDS data (A -I- B half-arrays). Also, in all 
cases the presence of the cut-off is if significant at greater 
than the 99.5% confidence level using the F-statistic. 

We have also tested whether the PDS so-called "spike" 
events can contribute to the PDS spectral shape. The an- 
swer is no: in our two observations of MR 2251-178 the 
contribution of the "spike" events to the PDS spectral 
shape is negligible. 

However, PDS event filtering using the so-called 
"Variable Rise Time" (VRT) selection does bring a sig- 
nificant improvement to the quality of fit, by increasing 
the signal to noise ratio of the PDS data. Therefore this 
filtering has been applied to all our PDS data. 

3.2.2. The pexriv model 

Our results with a pexriv model including a warm ab- 
sorber and an Fe Ka line (see also Fig. 0) are the following: 



a power-law slope F = 1.62 



+0.01 
-0.02' 



a high energy exponen- 



+38 
30 



tial cut-off of the incident power-law at Ecut-off = 115 
keV, the temperature of the ionized reflector Tdisk being 
fixed at 10^ K, an ionization parameter of the reflector 
£, = 465_24o ^^S cm s^ and finally the scaling factor for 



reflection Rion = 0.16lg'Jg, which is in fact the solid an- 
gle R=J7/27r subtended by the reflecting disk. It should 
be noted th at, for th ese reflection models, we use ^ as it 
is defined in R.FY99 , i.e. with the ionizing flux integrated 
between 0.01-100 keV. Because in the pexriv model ^ 
is defined between 5-20 keV we have made the neces- 
sary conversions, assuming a photon index F = 1.6 (i.e., 

CrFY99 ^ ^pexriv X 4.4). 

The fit statistics are good, with xt = 0-89 for 153 
dof. The low value of the scaling factor R tells us that 
most of the flux in this model is still coming from the 
direct, non-reflected component. Furthermore, one should 
bear in mind that a moderate ionization parameter such 
as obtained with this model, still produces a significant 
refiection "hump" . Since the present data obviously lack 
such a hump the scaling factor cannot be large. In all our 
reflection flts we flxed the inclination angle of the disk 
with respect to the observer so that cosz = 0.45. This was 
necessary because in our flts the viewing angle and the 
reflection scaling factor are not independent, nor are they 
well constrained with respect to one another at low values 
(R;^ 0.3 or cosi ^ 0.2). In fact, in the present pexriv flt 
the entire range of numerically permitted cos i values, i.e. 
cosz = 0.05 — 0.95 lies within the 99% confidence interval 
for this parameter. At cosi ^ 0.05 R is poorly constrained 
with a lower limit value of R> 0.4. At cos i ^ 0.95 we have 
0.05 < R < 0.22. R and cosi are linked in the fitting pro- 
cess because smaller values of i produce stronger reflection 
humps and vice-versa. 



3.2.3. The pexrav model 

Reflection on cold matter (e.g. neutral Gompton reflec- 
tion model pexrav in XSPEC) gives, compared to pexriv, 
a slightly worse fit with xt — 0.93 for 154. The difference 
is significant at the 99% confidence level. The neutral re- 
fiection scaling factor is weakly constrained with an upper 
limit Rneut. < 0.38. 



3.2.4. A self-consistent ionized reflection model: 
RFY99 

The pexriv model in XSPEC does not include Fe Ka emis- 
sion and one therefore needs to add a separate Gaussian 
line component to the fit model. Because of this draw- 
back we have also tested the self-consistent ionized reflec- 



tion models published by Ross and collaborators ( |RFY99| ) 
which include Fe Ka line emission and take into account 
the destruction of Ka photons by Auger effects. The mod- 
els were specially computed for the values of the incident 
power-law slope observed in MR 2251-178. As with pexriv 
we set the iron abundance at the solar value and included 
a warm absorber. The RFY models do not include the "in- 
cident" continuum so it is necessary to add it separately, in 
this case an exponentially cut-off power-law. We hereafter 
refer to this combination of components as the "RFY99 
model" . The power-law slope of the incident continuum 
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and the slope specifying the reflection component were 
fixed to have the same value. It should be noted that, un- 
like pexriv, the RFY99 reflection assumes a sharp, step- 
like cut-off 1^ of the incident spectrum at exactly 100 keV. 
This causes the RFY99 reflection component to decline 
steeply above ~50 keV. The RFY99 model gives a good a 
fit with xl — 0-94 for 155 dof , see Fig.|[ If the incident 
spectrum is taken to be a power-law without an exponen- 
tial cut-off the fit becomes significantly worse {xt — 1-40). 
The "RFY99" reflected fraction is ~0.11 and is compatible 
with the pexriv value. Both models also give comparable 
high energy cut-offs at ~100 keV. We have checked that 
this value is not an artefact caused by the sharp cut-off as- 
sumed in the RFY99 reflection component (at Eicst = 100 
keV and Eobs — 94 keV). Indeed, if a "hard tail" is added 
to the RFY99 reflection component, with an exponential 
cut-off at 200 keV, the fitted continuum cut-off energy 
remains at ~100 keV with a 90% confidence upper limit 
Ecut <200 keV. The reason for this is most likely the rel- 
atively low reflection scaling factor, R'-~^0.11. 

A noteworthy analogy can be seen between the pexriv 



and RFY99 models concerning the ionization parameter 
of the reflector. Due to non-linear x^ space two distinct 
flts are possible with the RFY99 models. In the first case 
the ionization parameter is ^ = 1625lg3Q erg cm s~^ and 
in the second ^ < 35.4 erg cm ■sT^ . This second value cor- 
responds to a nearly neutral reflector and gives a worse 
flt than with the first value (Ax^ = 4.4). Such a behav- 
ior is analogous to the pexriv versus pexrav solutions 
and therefore gives support to the validity of the ionized 
reflection models and their fitted parameter values. 

These reasons lead us to consider hereafter only the 
flt solution yielding the higher value of ^. The ionization 
parameter measured with pexriv is poorly constrained 
making it difficult to compare with the RFY99 value. 

The value of i (with the RFY99 model + 3 edges) 
is driven mainly by the two following factors: flrst, the 
intensity of the Fe Ka line and second, the low energy 
(<6.4 keV) spectral shape. This is verifled by performing 
the fits over limited energy ranges. For instance, excluding 
the 5.5-8 or 5-12 keV range leads to significantly higher 
ionization parameters, and weaker Fe lines. Likewise, if 
the range below 5 keV is excluded, then ^ decreases sig- 
nificantly. 

^ = 1625 implies a hot Fe Ka line in the reflection 
component, with Brest ~ 6.8 keV and Eobs ~ 6.4 keV. This 
is "hotter" than the line energies found with the other, 
non-self-consistent, fit models in Table 1. However, from 
Fig. H we see that the "hot" line is still consistent with 
the MFCS data. Indeed, no large fit residuals appear at 
the Fe Ka fine with the RFY99 model. 

Nevertheless, an improvement in x^ can be brought to 
the global RFY99 fit if an additional Fe Ka line (Eiost — 






Energy (keV) 

Fig. 5. Fe Ka line profile: MFCS data to model ratio in 
observer's reference frame. The model is a cut-off power- 
law with a Gaussian emission line for Fe Ka, where the 
line normalization has been set to zero for the plot. The 
energy of neutral Fe Ka emission (6.4 keV) is indicated 
by the dashed line 



As mentionned in the RFY99 paper, extending the illumi- 
nating spectrum to higher energies would raise the correspond- 
ing portion of the emergent spectrum via Compton downscat- 
tering of higher energy photons. 



6.53I^;i^ keV, EW = bbZH eV) is included (xt = 0.9, 
Ax^ = 8.8 for A(dof)=2). But the model is then no longer 
self-consistent and, in order to "compensate" for the addi- 
tional line, the ionization parameter of the reflector tends 
to increase and the reflected fraction to decrease. It may 
be a hint that a slightly different Fe abundance is required. 
For instance, increasing the Fe abundance increases both 
the line flux and the depth of the Fe K-edge. However, 
any iron over-abundance must be relatively low. This is 
because the ionization parameter derived from the RFY99 
model indicates that the reflecting matter is in the regime 
where ions are already too highly ionized to permit Auger 
destruction of iron line photons. Therefore these photons 
escape from the re flector to produce the observed hot iron 
line (Fabian et al. ^000 ). Another possibility is that the 
line is due to neutral reflection from matter further out, 
e.g. the torus or a cold phase of the accretion disk. 



3.2.5. The Fe Ka line 

Fits of the Fe Ka line (see Fig. ||) show that it is weak 
but very signiflcant at >99.9% confidence level (the F-test 
value for the presence of a line in the cut-off power-law 
model of Table |l| is F = 15.9 for a difference of 2 dof). 
A narrow line fit gives Eiine(i'est) — 6.53lo'x2 keV, an 
equivalent width EW = 62^25 ^V (the continuum is here 
the best fit pexriv model). The line energy is consistent 
with ionization stages Fe Xiii-xxiii. A significantly worse 
fit is obtained if one forces the line energy to be 6.4 keV 
(>95 % confidence level). In the flt models of Table || the 
width of the narrow line has been frozen at cr = 0.1 keV 
because no there is no significant improvement in the fit 
statistics when leaving a free. 
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Nevertheless, we checked whether broad wings are 
present beneath the narrow core. We fitted the MECS con- 
tinuum spectrum with a power-law continuum and galac- 
tic absorption to which we added, in turn, a narrow line 
or a "narrow plus broad Gaussian line" . In the second 
case both lines were required to have the same centroid 
energy. There is no strong evidence for the presence of a 
broad Gaussian component of the Fe K line. Such a broad 
line is not required by the data, the confidence level for 
its inclusion being only '^80%, using the F-statistic. The 
resulting line width is Abroad — l-SS^o^yy keV. The equiva- 
lent widths of the core and the broad components are then 
53.0 and 116.0 eV, respectively. So, we cannot completely 
rule out some broadening of the line by e.g. Compton scat- 
tering. 

Likewise, we have performed fits with the XSPEC 
diskline model using the MECS data alone and com- 
pared them with the broad gaussian line fits. The fit qual- 
ity of either model is statistically equivalent. Therefore 
it is not possible to formally exclude the presence of a 
relativistically "blurred" line. The width of the peak line 
emission remains lower than the instrument's energy res- 
olution at 6 keV (~8% in the MECS). 

When testing the diskline model we tried to min- 
imize the amount of free parameters by estimating the 
outer radius of the disk. We considered the simple case of 
an annulus illuminated from one side and a non-rotating 
central black hole with a mass of 2.5 x 10® Mq. This m ass 
was derived for MR 2251-178 by Brunner et al. ( |l997D by 
fitting geometrically thin a accretion disc models. They 
also found an inclination of the disc of 70 (39-82) degrees. 
From the value of the ionizing X-ray flux and the ionisa- 
tion parameter obtained with the RFY99 fit it is possible 
to estimate the surface area of the ionized disk, assuming 
a hydrogen number density of ng ~ 10^^ cm"'^. From the 
value of this surface area one can derive the outer radius 
of the disc (~6.46 Rg), the inner radius being fixed at 
6 Rg- The resulting diskline feature is double peaked, 
several keV broad and skewed. The blue peak is, however, 
about 4 times stronger in intensity than the red peak. And 
the FWHM of the blue peak is <0.1 keV (not resolved by 
BeppoSAX). 

In conclusion, our results with a narrow emission line 
confirm the previous measurements of the line made at 



higher energy resolution with ASCA (Euno = 6.57too? 
keV, EW = 62 ± 32 eV, a < 0.22 keV, Reeves et al. 1997) 
and they also are consistent with the EW-luminosity re- 
lation suggested by Iwasawa & Taniguchi (1993) and con- 
firmed by Nandra et al. ( |l997h| .) This relation is a strong 
anticor relation between the luminosity and the EW, an 
"X-ray Baldwin effect" which suggests Fe Ka line core 
EWs lower than ~80 eV at L2_iokoV ^ lO**^ erg s~^, 
which is the luminosity of MR 2251-178. 



3.2.6. The complex low-energy X-ray spectrum 

A soft excess component is not required by the data. The 
ionized gas of the warm absorber can account for the spec- 
tral complexity below 2 keV. Likewise, we see that the best 
fitting models (which include a warm absorber) do not re- 
quire excess cold absorption. The upper limit on excess 
neutral absorption is 0.16 x 10^" cm^^, using the RFY99 
model. 

The soft X-ray spectrum of MR 2251-178 is complex 
and can be well fit by considering the effects of strong 
absorption by partly ionized matter. Below 1 keV this 
absorption is essentially due to highly ionized oxygen. 
BeppoSAX cannot spectrally resolve the two oxygen edges 
O VII and O viii at 0.74 and 0.87 keV in the source rest 
frame. Therefore fits have been performed by multiply- 
ing the continuum model by either a single "blended" 
absorption edge or two oxygen edges with energies fixed 
at the values above. The "blended" edge parameters are: 
Ecdgc = 0.81 ± 0.01 keV, r = 0.80 ± 0.05. In the sec- 
ond case, T07 = 0.38 ± 0.07 and tqs = 0.48 ± 0.06. A 
third edge brings further significant improvement to the 
fit statistics (F-statistic F=21.41 for Adof=2). This third 
edge (E - 1.32 keV, r - 0.20) can be due to Ne x K, Mg 
i-iv K or Ni xiv-xv L. 

Alternatively, fits with warm absorber models were 
tested. The models were calculated using ION98, th e 199 8 
version of the photo- ionization code ION (Netzer 1996 ). 
They were combined with a cut-off power-law continuum 
and a gaussian Fe Ka line. The best fit gives a column 



density of the warm absorber Nw — (8.1_ 



1.2, 



xlO^ 



5+0.004 



and an oxygen ionization parameter Uoxy — 0.018_q qq3 
The oxygen ionization parameter Uoxy is the unitless ra- 
tio of ionizing photon flux to the electron density and is 
defined by integrating the rate of photon emission over 
the oxygen K-shell continuum, i.e. between 0.538-10 keV. 
In the present case Uoxy translates to an X-ray ioniza- 
tion parameter (defined by integrating between 0.1-10 
keV) Uo.i-io = Omotoml- The fit statistics with the 
warm absorber model are significantly better (at more 
than 90% confidence) than with the triple absorption edge 
fits: Ax^ = 5.13 for Adof=2 , leading to an F-statistic F 
= 2.9. 

We note that the ionization parameters of the reflec- 
tion component and the warm absorber are similar: using 
the apropriate conversion and flux integration between 
^0.01- 100 and Uq.i-io (e.g. Netzer 1996, George et al. 
|1998|) the reflector has Uo.i-io = 0.04(0.02 - 0.22), with 
the pexriv fit and Uo.i-io = 0.14(0.06 - 0.27), with the 
RFY99 fit, respectively. But the warm absorber is opti- 
cally thin and the RFY99 reflector optically thick. So, de- 
spite their similar ionization parameters they most likely 
are distinct. 

The comparison of our warm absorber fit results with 
previous warm absorber results for MR 2251-178 is com- 
plicated by the different formalisms used in the liter- 
ature to define the ionization parameter, in particular 
the integration range of the ionizing photon flux. For a 



8 
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start, we note that our measured warm absorber column, 
Nw — (8-lii'2) ^ 10^^ cm~^ is consistent with val ues de - 
rived from Ginga and EXOSAT data (Mineo et al. |l993|) , 



even when MR 2251-178 was in weaker or brighter states 
than seen by BeppoSAX. An earlier study of the same 
EXOSAT data by Pan et al. ( |1990[ ) seemed to suggest 
a column versus flux trend whereby the warm absorber 
column decreases with increasing X-ray flux. However, re- 
sults from Mineo et al. ( 1993 ) as well as from Otani et al. 
( |l997|) , for ASCA data, excluded such a trend. Unlike the 
EXOSAT data with relatively low statistics, the high sig- 
nal to noise ASCA spectra indicate slightly lower values of 
Nw than the BeppoSAX spectra do (< 7.5 x 10^^ cm"^). 
Finally, the Ginga, EXOSAT and ASCA data aU show a 
well defined correlation between the ionization parameter 
and the X-ray flux. Using the ionization parameter "trans- 
lation" curves of George et al. ( |1998| , Fig. 1) which give the 
relations between various definitions of the ionization pa- 
rameter, ^, U and Uo.i-io, we estimate that the ionization 
parameter of the warm absorber in MR 2251-178 globally 
agrees with the EXOSAT, Ginga and ASCA "U versus 
flux" correlation curves. To summarize, it seems that the 
warm absorber in MR 2251-178 is responding in a consis- 
tent way to the slow variations of the ionizing continuum 
flux and therefore is at most times in a state of ioniza- 
tion equilibrium. Any change in absorber column density 
between the 1996 ASCA and the 1999 BeppoSAX obser- 
vations could be explained by bulk motion of the absorber 
material across our line of sight. 



Whereas the value of the Fe Ka line energy, as mea- 
sured by ASCA and BeppoSAX in MR 2251-178, clearly 
excludes neutral Fe, the case for ionized reflection is less 
firm. The data strongly suggest that little reflection is 
present and that we are in fact mainly observing the direct 
ionizing continuum since very good fits are obtained with 
an exponentiall y cut-off power-law (x^ = 0.94, dof=:155) 
or a pexriv or RFY9£ model allowing the direct compo- 
nent to "shine" through at more than 80% (x^ = 0.89, 
dof=153 with pexriv, see also Table |l]). Because of the 
weak contribution of the reflection component the shape 
of the cut-off at ^100 keV must be intrinsic to the primary 
X-ray emission. However, despite its weak contribution, 
an ionized reflector can produce the observed ionized Fe 
emission, as shown with the RFY99 fit. 

Broad band observations of Syl galaxies with 
BeppoSAX have revealed exponential cut-off energies 
<,\QQ keV in only a couple Syls: e.g. NGC 5548 
(Ecut-off = 115^27 keV, Nicastro 

|1998D 



et al. pOQOD and 
Whereas in other 



NGC 4151 (70 keV, Piro et al 
Syls the cut-off occurs at much higher values, with lower 
limits around 200-300 keV (c.f. Gondek et al [l996t Matt 
1998). More recently, Zdziarski et al. (2000) have mea- 
sured the average OSSE spectra of Sy 1 and Sy 2. Using 
a cut-off power-law they derived expontenial cut-off en- 
ergies of 120l6o° and 130^50° keV for Sy 1 and Sy 2s, 
respectively. It appears therefore that the cut-off energy 
measured by BeppoSAX for the RQ QSO MR 2251-178 is 
compatible with the average OSSE Sy 1 value as well as 
the BeppoSAX values obtained so far for Sy Is. 



4. Discussion 

Both the XSPEC pexriv and the RFY99 models give very 
good fits to the BeppoSAX spectrum of MR 2251-178. Of 
all the models we tested only the RFY99 model can by 
itself and in a self-consistent way account for the ionized 
Fe line, part of the high energy cut-off and the lack of 
reflection hump, however, this model assumes a reflector 
with constant density. The upper layers of the reflector, 
which are subject to strong external illumination probably 
have lower densities and higher effective ionization param- 
eters than the inner layers. This would te nd to suppress 
Fe Kg em ission in the outer layers ( RFY99 ). Nayakshin et 
al. (200C) and Ballantyne et al. (2001) have recently cal- 
culated models of reflection on an ionized atmosphere in 
hydrostatic equilibrium. Nayakshin et al. (2000) flnd that 
in certain conditions Fe Ka photons can nevertheless be 
emitted: «/the incident spectrum is steep (F ^ 2) or if the 
intrinsic disk flux exceeds the X-ray illuminating flux, in 
these cases the reflecting layer is expe cted to yield ionized 
iron emission lines. Ballantyne et al. ( 200l| ) show that the 
constant density models of Ross & Fabian (1993) actually 
often give accurate flts to their "hydrostatic" reflection 
spectra, both statistically and parametrically, assuming 
an ASCA-like spectral resolution. On the other hand, the 
XSPEC pexriv and pexrav models give only poor fits to 
their spectra. 



5. Conclusion 

We have for the first time observed MR 2251-178 in the 
0.1-200 keV energy range using BeppoSAX. The spec- 
tral fitting described above shows that this source possi- 
bly combines two distinct signatures of ionized reflection 
which, to our knowledge, have not been detected together 
before in an individual radio quiet AGN: a strongly ion- 
ized Fe Ka emission line and a contribution, in this case 
relatively weak, from an ionized reflection continuum. A 
spectral break is apparent at higher energies which can 
be parametrized with the help of an e-folded power-law. 
The cut-off energy is measured to be in the range ~80- 
190 keV (90% confidence interval) and is compatible with 
values found for Sy 1 galaxies. 
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